Nonlinear generalized master equations and 
accounting for initial correlations 



Victor F. Los* 
Institute of Magnetism, Nat. Acad, of Sci. of Ukraine, 
36-b Vernadsky Blvd., Kiev 03142, Ukraine 

February 6, 2008 

Abstract 

By using a time-dependent operator converting a distribution func- 
tion (statistical operator) of a total system under consideration into the 
relevant form, new exact nonlinear generalized master equations (GMEs) 
are derived. The inhomogeneous nonlinear GME is a generalization of 
the linear Nakajima-Zwanzig GME and is suitable for obtaining both the 
linear and nonlinear equations describing the evolution of the (sub)system 
of interest. Actually, this equation provides an alternative to the BBGKY 
chain. To include initial correlations into consideration, this inhomoge- 
neous nonlinear GME has been converted into the homogenous form by 
the method suggested earlier in I10| . Obtained in such a way exact 
homogeneous nonlinear GME describes all stages of the (sub) system of 
interest evolution and influence of initial correlations at all stages of the 
evolution. The initial correlations are treated on the equal footing with 
collisions and are included in the modified memory kernel acting on the 
relevant part of a distribution function (statistical operator). No approxi- 
mation like the Bogoliubov principle of weakening of initial correlations or 
RPA has been used. In contrast to homogeneous linear GMEs obtained in 
|51llL)|. the homogeneous nonlinear GME is convenient for getting both a 
linear and nonlinear evolution equations. The obtained nonlinear GMEs 
have been tested on the space inhomogeneous dilute gas of classical par- 
ticles. Particularly, a new homogeneous nonlinear equation describing an 
evolution of a one-particle distribution function at all times and retain- 
ing initial correlations has been obtained in the linear in the gas density 
approximation. This equation is closed in the sense that all two-particle 
correlations (collisions), including initial ones, which contribute to dissi- 
pative and nondissipative characteristics of the nonideal gas are accounted 
for in the memory kernel. Connection of this equation at the kinetic stage 
of the evolution to the Vlasov-Landau and Boltzmann equations is dis- 
cussed. 
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1 Introduction 



Derivation of equations allowing for describing the evolution of the measurable 
(statistical expectation) values, which characterize a non-equilibrium state, re- 
mains one of the principal tasks of statistical physics. To get kinetic (or other ir- 
reversible) equations from underlying microscopic (reversible) classical or quan- 
tum dynamical equations, several approaches are usually used, which generally 
start with the Liouville-von-Neumann equation for a distribution function (sta- 
tistical operator) of the many-body system under consideration. 

One of them leads to the chain of coupled first-order differential equations 
for many-particle distribution functions known as the BBGKY hierarchy pQ. In 
order to obtain from the BBGKY hierarchy a closed kinetic equation for a re- 
duced (e.g. a one-particle) distribution function, one needs to decouple the chain 
in some approximation using a small parameter such as an interaction between 
particles or a density of the system n. Actually, it may be done by disregard- 
ing all initial (existing at the initial time moment to) correlations. The most 
consistent approach to decoupling the BBGKY hierarchy was developed by N. 
Bogoliubov PP. He suggested the principle of weakening of initial correlations 
which implies that at a large enough time, t — to ^> t cor (t cor is the corre- 
lation time caused by inter-particle interaction), all initial correlations damp 
and, as a consequence of that, the time-dependence of multi-particle distribu- 
tion functions is completely defined by the time-dependence of a one-particle 
distribution function (Bogoliubov's ansatz). The latter assumes an existence of 
the time interval 

tcor <C t — t() <C t re l (1) 

(t re i is the relaxation time for a one-particle distribution function) and leads 
to an approximate (valid only on the pointed above large timescale) conversion 
of the inhomogeneous (including two-particle correlations) equation for a one- 
particle distribution function of the BBGKY hierarchy into a homogeneous one. 
By employing this approach, Bogoliubov was successful in deriving the kinetic 
(Markovian) equation for a one-particle distribution function, particularly the 
Boltzmann equation, which describes the time-evolution with the characteristic 
time t re i but is unable to account for the initial stage of the evolution, to ^ t ^ 

tcor ■ 

The Bogoliubov principle of weakening of initial correlations is introduced 
as a following boundary (initial) condition to the BBGKY chain (e.g. for the 
classical physics case) 

lim U(t,t ) [F N (xi, . . .,x N ,t ) - f r (xi, ■ ■ .,x N ,t )] = 0, 

t— to — >oc 

N 

f r (xi,...,X N) t) = Y[Fi(Xi,t), (2) 

i=l 

where Fn(x%, . . . , xn, to) is an iV-particle distribution function at t = to, Xi = 
(i"i,pi) is a coordinate of the i-particle in a phase space and U(t,to) is the 
operator defining the evolution of all coordinates Xi (back) in time (see (lllfl ). 
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The first condition (JIJ implies that at i — to ^ tcor all correlations between the 
particles vanish and the second inequality in t — to <C t re i allows to neglect 
(in the first approximation in n) the time-retardation and substitute in J3J) 
U (t, to)f r (to) with U (i, to)f r (t). Thus, the decoupling of the BBGKY chain is a 
nonlinear procedure J2Jl which results in the closed nonlinear (like Boltzmann's) 
equation for a one-particle distribution function F\{xi,t) and introduces an 
irreversibility into the Liouville-von-Neumann equation. 

It should be noted, however, that actually only the small-scale initial cor- 
relations with t cor <C t re i may be disregarded on the kinetic timescale. Trying 
to get the collision integral in the kinetic equation with account for collisions 
of three and more particles (the terms of the second and higher order in n) one 
should deal with the divergencies caused, particularly, by the large-scale initial 
correlations with t cor > t re i (see, e.g. [2] and Thus, the correct expansion of 
the collision integral in powers of n and accounting for kinetic fluctuations are 
impossible if one does not account for initial correlations with t cor > t re i and 
Icor ^ I {Icor = vt cor , I = vt re i, v is the characteristic (mean) particle's velocity 
and I is the particle's mean- free path) 0]. 

In another approach, leading to the so called generalized master equation 
(GME) [SlEj and an iV-particle distribution function (or statistical operator) 
is divided by means of time-independent projection operators P and Q into the 
relevant and irrelevant parts 

F N (xi, . . . ,X N ,t) = f r (xi, . . . ,X N ,t) + fi(xi, . . . ,X N ,t), 

f r (xi, . . . ,x N ,t) = PF N (xi, . . .,x N ,t),fi(xi, . . .,x N ,t) = QF N (xi, . . .,X N ,t), 
P + Q = l. (3) 

Note, that the relevant and irrelevant parts depend on coordinates of all N parti- 
cles in contrast to the reduced distribution functions (like F\(xi)). The relevant 
part, which is mainly of interest, is, as a rule, a vacuum (slowly changing) part 
of a distribution function (statistical operator), i.e. the part with no correlations 
(like f r (xx, . . . , Xpf, t) in Applying the projectors P and Q to the Liouville- 
von-Neumann equation (see ©), one arrives at the GME, which is the exact 
inhomogcneous non-Markovian equation for the relevant part of a distribution 
function with a source (irrelevant part) containing all many-particle correlations 
at the initial moment to. It should be underlined, that the time-independent 
(linear) projection operators P and Q, used in this approach, commute with 
the time-derivative operator in the Liouville-von-Neumann equation for the dis- 
tribution function (statistical operator) of the whole system. Therefore, the 
procedure (0 is a linear one and transforms the linear Liouville-von-Neumann 
equation into the linear GME. To exclude initial many-particle correlations (a 
source), Bogoliubov's principle of weakening of initial correlations or simply 
the factorizing initial conditions (random phase approximation (RPA) which is 
incorrect in principle [S]) are commonly employed, and it results in the approx- 
imate closed linear homogeneous GME for the relevant part of a distribution 
function (statistical operator). The latter is then used for obtaining a linear 
evolution equation for an interesting reduced distribution function. 
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The abovesaid means that the proper treating of initial correlations is a 
key to obtaining the closed kinetic equations. To include initial correlations 
into consideration, the method (based on the conventional projection opera- 
tor technique) turning the conventional linear inhomogeneous generalized mas- 
ter equations (GMEs) into the homogeneous form has been proposed recently 
pillU). By introducing an additional identity for the irrelevant initial condition 
term, it has become possible to express this term through the relevant part of 
a distribution function. Thus, the irrelevant initial condition (inhomogeneous) 
terms in time-convolution (non-Mar kovian) GME (TC-GME) 00 Ej and time- 
convolutionless GME (TCL-GME) ^2 EL caused by the correlations at an ini- 
tial (at t — to) state of the whole system, have been transferred to the memory 
kernels of GMEs governing the evolution of the relevant part of a distribu- 
tion function. That resulted in exact time-convolution and time-convolutionlcss 
homogeneous generalized master equation (TC-HGME and TCL-HGME) for 
the relevant part of a distribution function (statistical operator) of a many- 
particle system which account for dynamics of initial correlations through the 
modified memory kernels. These equations, which are obtained with no loss 
of information, describe the evolution (influenced by initial correlations) of the 
relevant part of a distribution function at all timescales including the initial one, 
t — t < t cor , when the transition from the reversible evolution to the irreversible 
kinetic behavior (caused by the stochastic instability of the system's dynamics) 
is expected to occur. 

The obtained linear time-convolution homogeneous generalized master equa- 
tion (TC-HGME) has been tested for the space-homogeneous dilute gas of clas- 
sical 32 and quantum jTjO] particles. However, in order to get from the actually 
linear TC-HGME the desired nonlinear equations (particularly the Boltzmann 
equation) describing the evolution of classical or quantum gas of particles, the 
additional approximation disregarding the change in time of a one-particle dis- 
tribution function (statistical operator) during collisions should be performed. 
This is because the time-independent projection operator (for, e.g. the classical 
particles) [5j 

N 

P = ^\Fi{xiM))y^=i [dx 2 ...[ dx N (4) 

i=2 ^ ^ 

(V is the volume of the system) leads to the relevant part of an TV-particle 
distribution function 

N 

f r (xi, ...,x N ,t) = [J^Fi(x i ,to)]i<i(a;i,i), (5) 

i=2 

which is linear with respect to the time-dependent one-particle distribution func- 
tion of interest Fi(xi,t) (Fn satisfies the normalization condition 

J dx\ . .. J dxjv-Fjv^i, ■ ■ ■ ,XN,t) = 1). Then, in order to get from ob- 
tained linear equation for Fi(x\,t), the nonlinear one, F\{xi,to) may be ap- 
proximated by Fi(xi,t) if t — t <C t re i (see QJ). Thus, in this approach we 
are not restricted by the first inequality l(Tjl. t cor <C t — to, but still need its 
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second part. It is difficult to imagine how to get a nonlinear symmetric in 
time relevant distribution function f r (xi, . . . , xjv,i) © by applying the time- 
independent projection operator P to the iV-particle distribution function. 

In this paper we will obtain a new class of evolution equations for the rel- 
evant part of a distribution function (statistical operator) suitable for deriving 
both the nonlinear and linear evolution equations for the reduced distribution 
functions (statistical operators) of interest and which are valid at all timescales. 
It will be exact nonlinear inhomogeneous and homogeneous generalized mas- 
ter equations deduced from the Liouville-von-Neumann equation with the help 
of (generally) nonlinear time-dependent operator P(t) converting the many- 
particle distribution function (statistical operator) into the appropriate relevant 
form. The time-dependent operator P(t) is not generally a projection operator. 
For the case of a time-independent linear projection operator P the obtained 
equations reduce to the linear conventional GME and obtained recently HGME. 
In contrast to these linear equations, the obtained nonlinear equations are more 
general and, therefore, more convenient, particularly, for studying the space- 
inhomogeneous nonideal gas of particles when the space-nonlocality and time- 
retardation are essential (see, e.g. |2|). In addition, the initial correlations may 
be exactly included into consideration by converting inhomogeneous nonlinear 
GME into the homogeneous form. 

At first (Sec. 2), we will derive the exact inhomogeneous nonlinear time- 
convolution GME for the relevant part of a distribution function (statistical 
operator) containing an irrelevant part (comprising initial correlations) which 
poses a problem do deal with as in the case of conventional linear GME. Ac- 
tually, this equation provides an alternative to the BBGKY chain. Like linear 
inhomogeneous GMEs, the obtained nonlinear equation will be transformed in 
Sec. 4 into the homogeneous nonlinear time-convolution GME by the method 
suggested earlier in [5]. The derived exact homogeneous nonlinear equation for 
the relevant part of a distribution function (statistical operator) accounts for 
initial correlations by means of the modified memory kernel of inhomogeneous 
nonlinear equation. It is worth noting that there is no way yet to exactly account 
for initial correlations dealing with the BBGKY chain. 

The obtained nonlinear GMEs will be applied to the space inhomogeneous 
nonideal dilute gas of classical particles by choosing the appropriate nonlinear 
operator P(t). In the linear approximation in the small density parameter (only 
two-particle correlations (collisions) are taken into consideration), the inhomo- 
geneous nonlinear GME results in the inhomogeneous nonlinear equation for 
a one-particle distribution function containing the self-consistent Vlasov term, 
the collision integral accounting for all space- and time-nonlocalities and an ir- 
relevant part (a two-particle correlation function at the initial moment of time 
to). Accepting RPA for t — to (or Bogoliubov's principle of weakening of initial 
correlations for t — > — oo) and disregarding the terms accounting for the space- 
and time-change of a distribution function in the collision integral on the mi- 
croscopic scales l cor and t cor (which are of the small density parameter order), 
we obtain the Boltzmann equation (Sec. 3). 

Then, in order to include initial correlations into consideration, we apply the 
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derived homogeneous nonlinear GME to the same system (Sec. 5). In the linear 
approximation in the density parameter we obtain the homogenous nonlinear 
equation for a one-particle distribution function valid at all timescales (no re- 
strictions like Q are needed) and accounting for influence of the dynamics of 
initial correlations at all stages of the evolution process. In this linear (in n) 
approximation, the evolution of collisions and initial correlations is governed by 
the exact two-particle propagator, and in this sense the obtained equation is 
closed. If all initial correlations vanish at the kinetic timescale (and all non- 
localities of the collision integral contributing to the thermodynamic functions 
of the nonideal gas are omitted) , this equation switches to the nonlinear Boltz- 
mann equation with the self-consistent nonlinear Vlasov term. The results are 
summarized in Sec. 6. 



2 Inhomogeneous nonlinear generalized master 
equation 

As it was already mentioned, the aim of this paper is to generalize the method 
suggested in jH] and jTH] for deriving the homogeneous generalized master equa- 
tions (HGMEs) retaining initial correlations to the case when the (projection) 
operator selecting the relevant part of a distribution function (statistical oper- 
ator) depends on time. It will allow for obtaining more general equations (in- 
homogeneous and homogeneous nonlinear generalized master equations) which, 
particularly, are more flexible in selecting the suitable for the problem under 
consideration relevant distribution function (statistical operator) describing the 
evolution of a (sub)system of interest. We will also generalize our approach to 
the case of the time-dependent external forces. 

We start with the Liouville-von-Neumann equation describing the evolution 
of a distribution function (statistical operator) F(t) of a total system under 
consideration 

^F(t)=L(t)F(t), (6) 

where, L(t) is the Liouville (evolution) operator acting on F(t) and depending, 
generally, on time. For example, in the case, when F(t) is related to the system 
of N classical particles and depends on their coordinates and momenta Xi — 
(fi,pi) (i=l,...,N), 



L(t)F(t) = {H(t),F(t)} p =J2 



A r 8H(t)dF(t) 8H{t)dF{t) 



dri dpi dpi dv t 



(7) 



where H(i) is the corresponding time-dependent Hamilton function of all X4. 
For the quantum physics case, L(t) is the superoperator acting on a statistical 
operator as 

L(t)F(t) = ±[H(t),F(t)}, (8) 
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where [,] is a commutator and H(t) is the energy operator (Hamiltonian). The 
distribution function (or statistical operator) is subjected to the normalization 
condition 

dx\ . . . I dxjvF(xx, . . . acjV) t) = 1 (9) 



or 



TrF(t) = 1. (10) 
The formal solution to equation may be written as 



F(t) = U(t, t )F(t ), U(t, t ) = Texp 



JdsL(s) 

to 



(11) 



where T denotes the chronological time-ordering operator which orders the prod- 
uct of time-dependent operators such that their time-arguments increase from 
right to left and F(to) is the value of a distribution function (statistical operator) 
at some initial time to. 

Let us introduce the generally time-dependent operator P(t) converting the 
distribution function F(t) into a relevant distribution function f r (t), which, as 
a rule, is a vacuum (with no correlations) many-particle distribution function, 
i.e. a product of the reduced distribution functions which are necessary and 
sufficient for describing the evolution of the measurable (statistical expectation) 
values. 

However, the problem is to obtain the closed evolution equation for such 
reduced distribution functions in view of the existence of correlations between 
the elements of the considered system including initial correlations at the initial 
time to. 

Thus, we define the relevant distribution function (statistical operator) as 

f r (t) = P(t)F(t) = P(t)U(t, t )F(t ). (12) 

Note, that generally the operator P(t) is not a projection operator. Moreover, 
the action of P{t) on F(t) is not a linear operation anymore (in contrast to 
applying the conventional time-independent projection operator [H]) because 
P(t) may depend on the time-dependent distribution function of interest (see 
below). Therefore, by applying P(t) to the linear Liouvillc equation ©, one will 
generally obtain a nonlinear equation. Because P(t )F(t) ^ f r {t) (at t ^ t), 
operator P(t) does not commute with the derivative (which is the case for 
the linear operators). 

Thus, using the Liouville-von-Neumann equation |j()J|, the equation of motion 
for f r {t) H12J) may be written as 



dfrjt) 

dt 



am ■ mm 



[fr(t)+fi(t)]. (13) 



Here, fi(t) = [F(t) — f r (t)], dP Jf > denotes the operator obtained by taking the 



derivative of the operator P(t), i.e., we assume that the operator dP Jf > exists. 
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To make sense of dP Jf > , one may define P(t) as 



P(t) = C(t)D, 



(14) 



where C(t) is a well defined (operator) function of time and D is a time- 
independent (super)operator acting on F(t) and integrating off all unnecessary 
variables in the spirit of the reduced description method (more details on the 
reasonable specification of P(t) will be given below). 

As it follows from (|13(l . the equation for the irrelevant part of a distribution 
function /j(t) = F(t) — f r (t) may be obtained by using equations JHJ and IJ13J1 
and splitting F(t) into the relevant and irrelevant parts. Thus, 



dfijt) 
dt 



Q(t)L(t) 



dP(t) 



dt 



[fi(t)+fr(t)], 



(15) 



where Q(t) = 1 - P(t) and, therefore, /;(*) = Q(t)F(t). 
The formal solution to equation l|15|) is 



fi(t) = fdtS(t,t) 

to 



Q(t )L(t ) 



dPjt) 
dt' 



where 



S(t, to) — T exp < Jds 



Q(s)L(s) 



fr(t) + S(t,t )fi(t ), 

dP(s) 



ds 



(16) 



(17) 



The latter operator may be presented as a following series 



S(t,t ) = 1 + Jdh 

to 

t t 
+ Jdt x J dt 2 

to t\ 



Q(h)L(t-_ 



Q(*i)i(*i) 
Q(t 2 )L(t 2 ) 

. dP(h) 



dP(h) 



dP(t 2 ) 



dt 2 



dh 



(18) 



Note, that if the operators P and L do not depend on time, then S(t, to) 

exp[QL(t-ta)]. 

Inserting (|TB|) into (|T^|l . we get 



dfrjt) 

dt 



P(t)L(t) 



dP(t) 



+ JdtS(t,t) 



dt 

Q(t)L(t 



{fr(t) 



dPjt 
dt' 



fr{t) 



S(t,t )fi(t )}. 



(19) 



Equation 119|) is the generalization of the Nakajima-Zwanzig time-convolution 
generalized master equation (TC-GME) (see jSlEllZj) for the relevant part of a 
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distribution function (statistical operator) to the case of time-dependent oper- 
ators P and Q. As it has been already noted, P(t) and Q(t) are not generally 
the projection operators in the usual sense and in our derivation of l|19|) we have 
never used such a property of the projection operators as P 2 = P, Q 2 — Q (for 
P(t)P(t ) and Q(t)Q(t ) such a property holds only for t = t , as it will be seen 
below). The obtained equation, like conventional TC-GME, is an exact inho- 
mogeneous integra-differential equation for the relevant part of a distribution 
function (statistical operator) containing the irrelevant part of a distribution 
function (a source) at initial moment of time /j(to) = F(to) — f r (to). But, 
equation l|19fl is generally the nonlinear equation for a reduced distribution 
function of interest (e.g. for a one-particle distribution function considered in 
the next section) and, therefore, is also convenient for studying the evolution 
of many-particle systems, which are described by the nonlinear Boltzmann-type 
equations (in the kinetic regime) and where the nonlinearity is caused by the 
collisions between particles. 

For time-independent P and Q equation (|19|) reduces to the conventional 
linear TC-GME. This type of equations is more suitable for studying, e.g. the 
evolution of a subsystem interacting with a large system in the thermal equilib- 
rium state (a thermostat). 

To simplify equation (|19|l and make sense of the derivative dP Jf > , let us 
consider the operator P{t) of the form (|14|l . P(t) = C(t)D. This representation 
of P(t) is suitable for studying many cases of interest (see the next section). 
Also, splitting a distribution function (statistical operator) into the relevant 
P(t)F(t) and irrelevant Q(t)F(t) parts makes sense if the operator P(t) satisfies 
the following relation 

P(t)P(t) = P(t) (20) 
for any t. Given the form of P(i), P(t) = C(t)D, equation (|2U[1 entails 

DC(t) = 1 (21) 

for any time moment t. Equation (|21|l generally represents the normalization 
condition for the distribution function(s) making up C(t) (see |J5J, <|10[) and 
H38[) below) and implies that C(t) depends on the variables which are removed 
from F(t) by operator D (f r (t) — P(t)F(t) depends on the complete set of the 
variables of a distribution function F(t) as indicated by 10). From the property 
l|21|) we have 

P(t)P(t) = P(t),Q(t)Q(t) = Q(t'),P(t)Q(t) = 0, 

Q(t')P(t) = P(t)~P(t). (22) 

The condition (|21|l also leads to the following relations for the time-derivative 
ofP(t) 

dP(t) dC(t) , . 

where we assumed that the operators D and d/dt commute (D does not depend 
on time). 
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Let us also introduce the reduced distribution function (statistical operator) 

fred(t) 

f red (t) = DF{t) = Df r (t), (24) 

which is actually needed for calculation of the expectation values of interest in 
the non-equilibrium state. Also, due to (|21[) . 

D dfr& = dfr^t) Dm = DQ{t)F{t) = Q_ (25) 

Thus, by applying the operator D to equation (|19fl from the left, we get 



<U "" ,:) =DL(t){f r (t) + Jdt'S(t,t) 



9t 



Q(t)L(t) dP ^ 



fr(t) 



Of 

S(t,to)fi(to)}, (26) 



where (|2*T|) . and (|23|l have been used. 

Without loss of generality the Hamiltonian of a system may be split into 
two terms as H(t) = H°{t) + H (t) (L{t) = L°(t) + L{t), correspondingly), 
where H°(L°) is related to the energy of the noninteracting particles (or sub- 
systems) and H (L ) describes the interaction between particles (subsystems). 
Then, except satisfying the condition J2U, the operator D, selecting the reduced 
distribution function, should commute with L°in the following manner 

DL°(t) - L° ed (t)D, (27) 

where L° ed (t) is the reduced Liouvillian L° depending only on the variables of 
the reduced distribution function (statistical operator) f re d{t), and l(77|) means 
that both sides of this relation act on the arbitrary function (operator) de- 
fined on the phase (Hilbert) space. Taking into account that C(t )L° ed (t) = 
L° ed (t)C(t ) (because C(t ) depends on the set of variables that are removed 
from F(t) by D and L° ed (t) does not depend on them), equation l|27() may be 
written in the more general form 

P(t)L°(t) = L° red (t)P(t). (28) 

The commutation relation (|27|l follows from a self-consistency reason, because 
if L (t) = 0, the evolution of a "free" (noninteracting with the rest of a system) 
subsystem should be described in equation (|26|l only by L° ed (t) (DL°(t)f r (t) = 
L° ed (t)f red (t)) (see also Q]). That this is the case may be seen from the relations 

DL°(t)S(t,t 1 )Q(t 2 ) = 0, 

DL°S(t,t 1 )^l=Q, (29) 

which follow from (|2*7Jl , lfT7|l , and JjZHJl . Thus, all terms in are propor- 
tional to L (t) but the first one (L® ed (t)f red (t)). 
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In the same way (using l(27|) . H17|) and relations that follow from (12111 ) the 
following equations may be proved 



S&tJQfa) = U(t,t 1 )Qfo), 

S(( , (l) ^=t7 (( , (l ) a ^> 



di 2 



(30) 



where 



Z7(t, ti) = Texp<^ jds L°(s) + Q(s)L 



(31) 



Thus, making use of ((29(1 and l|3U|) in 1)26(1 . equation (|19|l may be finally rewrit- 
ten as the following equation for the relevant part of a distribution function 
(statistical operator) 



dDfrjt) 

dt 



= DL(t)f r (t) 

t 

+ DL it) f dt U(t, t) 



Q(t )L(t ) - 



dPjt ) 
dt' 



fr{t 



DL (t)U(t,to)Mto), 



(32) 



where Df r {t) = f re d(t) and holds. 

The fulfillment of the basic conditions <(21(1 and 1(27(1 will be demonstrated 
in the next section for a gas of interacting classical particles. 



3 Dilute inhomogeneous gas of classical parti- 
cles 

Let us apply equation 1(32(1 to the case of a gas of N (N ^> 1) identical classical 
particles. For such a system the Liouville operator L may be presented as 



L = L° + L 



JY 



«=i 

JV 



__ 

rn 



Vi = 



dxi 



i<j=l 



(33) 



N 

where L° corresponds to the kinetic energy H° = pf /2m of particles with the 

8=1 

momenta Pi and mass m, while L corresponds to interaction between particles 

/ N 

H = ^2 Vij with the pair potential Vij = V{ | x^— Xj |). Thus, in the considered 

i<j=l 

case L does not depend on time. We do not suppose here that the interaction is 
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weak and assume that all necessary requirements to the properties of forces, by 
which the particles interact with each other, are met (particularly, the bound 
states are not formed). 

We will also assume, as usual, that all functions $(xi, . . . , Xn), defined on the 
phase space, and their derivatives vanish at the boundaries of the configurational 
space and at pi = ±00. These boundary conditions and the explicit form of the 
Liouville operators l|33|) lead to the following relations 

J dx i L° i $>(x 1 ,...,x N ,t) = 0, 

J dxi J dxjL i j^(xi, . . . , xn, t) = 0. (34) 

We are looking for an evolution equation for a one-particle distribution func- 
tion 

F 1 (x i ,t) = vjdx 1 ...J dx^ J dx i+1 • • • / dx N F( Xl ,...,x N ,t). (35) 
To this end, it is convenient to define the operator P(t) as 

N 



P(t) = C(t)D,C(t) 



.i=2 



d = y N _-i I "'•■"2 ••• / <i.r,\ ■ 



where V is the volume of the system. 

Generally, the s-particle (s ^ N) distribution function is defined as 



F s (xi, ...,x s ,t) = V s J dx s+1 ... j dx N F( Xl , ...,x N ,t) (37) 
and satisfies the normalization condition (sec ©) 

dxi... I dx s F s (x u ...,x s ,t) = V s . (38) 



Because the distribution function (|37|l for s — N, -Fjv(i) = V F(t) (hereinafter 
the dependence on Xi is often suppressed for brevity) , satisfies the same Liouville 
equation 10, we may apply the projection operator 13tj|) to F^(t) (instead of 
F(t)) and obtain the following relevant part of the distribution function 

f r (t) = P(t)F N (t) = flFi(*i,t). (39) 

i=i 

From the definition (|39l) one can see an advantage of the time-dependent op- 
erator P(t): by applying it to an iV-particle distribution function Fjy(t), it 
is possible to define the relevant part H39|) of a distribution function, which 
looks more natural and symmetric than used in |S] and obtained with the 
help of the time-independent projection operator |@J (we had to use the time- 
independent projection operator because the TC-GME considered in j^] is valid 
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only for such operators). Moreover, the relevant part l)39|) has a clear physical 
meaning because it represents a slow changing vacuum (with no correlations) 
part of an TV-particle distribution function which is mostly of interest. 
Thus, the irrelevant part 

fi(t) = Q{t)F N (t) = F N (t) ~ UFx(xi,t) (40) 

i=l 

describes all inter-particle correlations, particularly the initial ones (existing at 
initial moment of time to), and makes up the inhomogeneous term of 132|) . The 
irrelevant part of the iV-particle distribution function at to, fi(to) — Fjy(to) — 
f r (to)j may be always represented (as well as at any other time t) by the following 
cluster expansion 

N N-2 N N-3 

/i(*o) = X! 92{x l ,x j ) Yl Fi(x k )+ ^2 93(Xi ) n *!(*»)+•», 

i<j — l k^ij i<j<k—l l^i,j,k 

(41) 

N-2 N-3 

where J | Fi(xk) and J \ Fi(xi) stand for the products of N — 2 and N — 3 

one-particle distribution functions with k ^ i,j and I ^ i,j,k, respectively, 
whereas g2(xi,Xj) and g-$(xi,Xj,Xk) are the irreducible two-particle and three- 
particle correlation functions (further terms in (|41|l are defined in the same 
way). Here and further, the one-particle distribution and correlation functions 
with the suppressed time argument denote the functions taken at t = to (e.g., 
Fi(x k ) = F 1 (x k ,t )). 

Using the normalization condition (|38|) for a one-particle distribution func- 
tion, J F\{xi, t)dxi = V, it is easy to see, that by applying the integral operator 
D to C(t) (both defined by (HP) we obtain DC It) = 1, i.e. the condition J2TJ 
and the consequences of this condition given by (1221) and (|23[) hold. Also, in 
this case 

fred{t)=Df r {t)=F l {x 1 ,t), 

DQ(t) = 0,Dfi(t)=0 (42) 
in accordance with l|24|l and 1)25(1. 

The commutation rules (O and J2H1) also follow from (jSHJ , JHH) an d CSJ) 

DL° = L\D 

P(t)L a = L°P(t), (43) 

i.e. in the case under consideration L® ed = L\. Therefore, the identities (|29|) . 
<|3U|) and l|31|) also hold, where LP and L do not depend on time in the considered 
case. 

Making use of (|3"4*|) . (|4*2*|> and equation (|3*2*|) may be presented in the 

case under consideration as 
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dt 



L?Fi(n,i)+n J dx 2 L 12 F 1 (x 2 ,t)F 1 (x 1 ,t) 



t 



DL J dt U{t, t ) 



'0 rDn^V 



(L u -L») + Q{t)L - 



or 



fr(t) 



to 

+ DL'u(t,t )fi(t ), (44) 

where n = N/V is the density of particles, D, f r (t) and fi(to) are defined by 
, (jSfll and l|4*T|l , respectively. The evolution operator (fTTTTl may be presented 
by the series like (|T%j) 

t 

U(t,t) = 1 + /dti[i° + Q(ti)i] 
t' 

t t 

+ fdt 1 fdt 2 [L + Q(t 2 )L'}[L + Q(t 1 )L'} + ... (45) 
t' «i 

Let us consider equation (|44|) in the first approximation in the particles' 
density n. The corresponding dimcnsionless small parameter of a perturbation 
expansion is 

7 = rgn<l, (46) 

where ro is the effective radius of an inter-particle interaction. From the second 
equation (|34() it is easy to see, that all terms in the expansion l|45|l containing 
P(t)L result in the expressions, proportional at least to the first power of n (like 
the second term in the right-hand side of 144|l h Thus, in the first approximation 
in n all terms of (|45|l with P(t)L may be disregarded, because the terms of 
equation (|44l) containing U(t, t ) are already of the first order in n (application of 
DL leads to the expression of at least the first order in n) . In this approximation 
all terms under the integrals in (|45|l do not depend on time and the operator 
U(t,t ) reduces to 

U(t,t)=exp[(L° + L)(t-t')]. (47) 

Now, it is convenient to use the following expansion for the exponential 
operator (|4*7jl 

t 

e (L°+L')(t-t ) = e L\t-t) + f dee L a (t-e) L ' e (L«+L')e^ (48) 



Making use of this expansion and relations (|34fl . we get from (|44f) the following 
evolution equation for a one-particle distribution function in the linear approx- 
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imation in the density parameter l|46|) 



dF 1 {x 1 ,t) 
dt 



= L° 1 F l (x 1 ,t) + n J dx 2 L 12 F 1 {x 2 ,i)F 1 (x 1 ,t) 

+ n [ dx 2 L' 12 t ] t °dt 1 e^ +L '^ tl [(L° 2 + L 12 )F 1 (x 2 ,t - t^F^t - t x ) 



o 




where L\ 2 = Li + L 2 . Obtaining (|49|l . we have also used the definition for the 
irrelevant part of a distribution function (|41|) and taken into account that each 
additional integration over 23, . . . adds an additional power of n (therefore, in 
the linear approximation in n all formulae contain no more than one integration 
over the phase space). 

Equation (|49|l is the main result of this section. It is necessary to stress that 
equation l|49[l is a nonlinear inhomogeneous time-convolution (non-Markovian) 
master equation for a one-particle distribution function containing initial cor- 
relations (a source). Such an equation cannot be derived from the conventional 
linear time-convolution Nakajima-Zwanzig generalized master equation (TC- 
GME) [SI El EI j which follows from the linear Liouville equation with the help of 
the time-independent linear projection operator. Equation (|49|l is exact in the 
linear approximation in the small density parameter 1(46(1 and is valid on any 
timescale and for any inhomogeneity of the system under consideration. The 
first term on the right-hand side is a conventional flow term. The second non- 
linear Vlasov term represents the self-consistent field acting on the given (first) 
particle and defined by all particles of the system. It is worth underlying an 
impossibility of getting (without additional approximation) this term (showing 
up only for a space inhomogeneous case) from the GME (or HGME) containing 
the time- independent projection operator (see [5]). The term, quadratic in L 12 , 
describes the collisions between particles and leads to dissipation in the system. 
The terms with L 2 (derivative in space) and with the time-derivative account 
for the change in space and time of the one-particle distribution function on the 
microscopic (see also below) scales of the orders of r$ and t cor ~ ro/u, respec- 
tively (v is the mean particle velocity). These two terms define the contribution 
of the pair collisions to the nondissipativc characteristics (thermodynamic func- 
tions) of the non- ideal gas (see, e.g. 0]) and are usually absent in the standard 
derivation of the kinetic equations in the linear in n approximation. These terms 
usually appear in the next (second) approximation in n (see, e.g. |13p that is 
not very consistent with the dissipative three-particle collisions term showing 
up in this approximation. The last (irrelevant) term in l|49|l accounts for initial 
correlations (at t = to) and in the linear approximation in n is given by the 
two-particle correlation function (72(2:1, x 2 ). 

Let us look at the terms of 149(1 more closely, especially at their behavior 
in time. First of all, evolution with time is governed by the exact two-particle 
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propagator G\2(t) — exp[(L°2 + L 12 )t] (it is natural for the considered dilute gas 
in the lowest approximation in density), which satisfies the integral equation 

t 

G 12 (t) = G° 12 (t) + y#i<^ 2 (t-ti)L' 12 Gia(ti), (50) 
o 

where G\ 2 {t) — ex P(^i 2 )t is the propagator for non-interacting particles. The 
action of G\ 2 (t) on any function defined on the phase space is given by 



e L *$(xi, ...,x N ,t) 



N 



n 



$(xi, . . .,X N ,t) 



e 

.i=i 

= $(xi - Vit, pi, . . . , xjv - VNt, Pat, t), (51) 



which follows from (|33|) . Thus, if the particles' dynamics possesses the necessary 
properties such as, e.g. a mixing ergodic flow in a phase space, caused by the 
local (stochastic) instability, it may be expected that under the action of Gn (t) 
the distances between particles will rapidly increase with time (a contribution of 
a "parallel motion" is negligible). If the effective interaction between particles 
V(\ri — vanishes at a distance |r< — Tj\ > ro, then, the integrand of the 
integral over t\ in l|49l) vanishes at t — to > t cor and the same is generally (but 
not always) valid for correlation function 52(^1, x%) which also depends on the 
interaction between particles (in more detail it has been considered in 9 for 
the space homogeneous case). 

On the other hand, the condition 146|) implies the existence of the time 
hierarchy 

tcor ^ trel: (52) 

where t re i ~ r )~ 1 t cor is the relaxation time for a one-particle distribution func- 
tion Fi(xi,t) (1/Wq is a particle's mean free path). In the initial stage of the 
evolution f t ^ t cor , which is very interesting and essential for understand- 
ing the irreversibility problem and for studying the non-Markovian processes 
(the memory effects) and the ultrafast relaxation effects, the initial correlations 
may be important. However, if one considers the timescale t — to 3> t cor , the 
initial correlations may damp and the source in (|49|l containing g%{xx,X2) may 
be disregarded. But as Bogoliubov noted [Q, this can be only done if the time 
interval £Q) exists due to 1)52(1 . The latter conclusion follows from the fact that 
the corrections to the solution of 149fl . due to initial correlations, are propor- 
tional to \t — to | (secular terms) and are small only at \t — to\ <§C t re i (see 0). 
Thus, the expansion of the irrelevant initial condition term in n (as we did it 
while going from (|44|l to 1(49(1 ) is ineffective because the terms of this expansion 
are proportional to time and the convergence is very poor at the timescale of 
interest in the kinetic theory t — to > t re i . 

To get rid of the undesirable initial condition terms, it is often assumed that 
fi(to) = 0, which corresponds to the random phase approximation (RPA) and 
is incorrect in principle (see |Hj). 
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Therefore (as it is already mentioned in the Introduction), in order to study 
the kinetic stage of the evolution at t > t re i a more sophisticated approach like 
that based on the principle of weakening of initial correlations (Bogoliubov's 
anzats) ,1. is needed. However, the Bogoliubov approach does not allow for 
considering the initial stage of the evolution to ^ t ^ t cor and the situations 
when the large-scale correlations, associated, say, with the quantities that are 
conserved, may exist. These correlations may only vanish on the timescale of the 
order of the relaxation time t re \ . This approach also implies that the correlations 
with t cor > t re i appearing in the higher approximations in n are not essential, 
although they are essential (see, e.g. 00]). 

If F\ (xi , t) changes significantly on the macroscopic (hydrodynamic) length 
lh and the corresponding timescale th = lh/v (for the space homogeneous gas 
there are only the characteristic scales I and t re i on which F\ (xi , t) changes) the 

terms with L\ and dFl ( x ^~ l ) m equation l|49|) may be disregarded, because 
they are both at least of the order of 7, i.e. their contribution to dissipative 
processes is of the second order in n (however, these terms should not be omitted 
if we wish to account self-consistently for the nondissipative characteristics of 
the nonideal gas). At the kinetic timescale, t — to > t re i, the upper limit of 
integration over t\ may be extended to infinity (to — * —00). If we also accept 
the RPA for to — > —00, i.e. put 32(^1,2:2) = 0, then equation (|49|l reduces to 



where we have taken in consideration that Fi(xi,t — t\) may be substituted 
(within the adopted approximation) with F\ (x^ ,t) at t > t re i 3> t cor (the in- 
tegration over t\ in 149[) gives essential contribution only up to t\ ~ t cor ). In 
the last (collision) term of l|53|) the distribution functions of the collided parti- 
cles may be taken at the same space point ri = V2 with the adopted accuracy 
~ ro/r h < 1 (or r Q /l < 1). 



If the integral J exists (and this is stipulated by the properties of the par- 



ticles dynamics), equation l|53|l represents the nonlinear Markovian irreversible 
in time kinetic equation. In the case of weak interaction between particles this 
equation in the first approximation in the small interaction parameter e gives 
the reversible Vlasov equation (equation l|53|) with disregarded the third collision 
term in the rhs, which is of the second order in e). In the second approxima- 
tion in e it coincides with the Vlasov-Landau kinetic equation (see, e.g. |14|1. 
For the space homogeneous case, when a one-particle distribution function does 
not depend on a particle coordinate, F 1 (xj,t) = Fi(pj,t) (J Fi(p,t)dp —1), the 
second (Vlasov) term in (I53|) vanishes, because we consider a potential Vij lKS.it 
dependent on the particles' coordinates difference. In this case equation l|53|) is 
equivalent to the Boltzmann equation (see, e.g. Jl]). It should be underlined, 
that the nonlinearity of l|53|l appears quite naturally in this approach (based on 



dF 1 {x u t) 
dt 




(53) 



00 
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time-dependent P(t) ffify ) in contrast to the time-independent projection oper- 
ator approach (see [HIE])- In the latter case, we had to approximate Fi(x2,to) 
by Fi(x 2 ,t) (it is possible on the timescale given by Q) in the evolution equa- 
tion in order to get a nonlinear equation from the linear GME or HGME. 

4 Homogeneous nonlinear generalized master equa- 
tion retaining initial correlations 

In section 2, the nonlinear generalized master equation for the relevant part of a 
distribution function has been obtained. However, equation l|19|l is not closed in 
a sense that it is inhomogeneous one and contains an irrelevant term comprising 
generally all multi-particle correlations. As was pointed above, to get rid of this 
irrelevant part and obtain the closed equation for, e.g., one-particle distribution 
function (like JSSJl), the Bogoliubov principle of weakening of initial correlations 
or RPA are usually used. 

In order to include initial correlations into consideration and to get the evo- 
lution equations, valid on all timescales, the method turning the conventional 
linear inhomogeneous generalized master equations (GMEs) into the homoge- 
neous form has been recently proposed in This method has been applied 
to conventional linear GMEs (obtained by means of the time-independent pro- 
jection operator technique) and that resulted in the linear homogeneous general- 
ized master equations (HGMEs). Naturally, these HGMEs are more convenient 
for studying the linear situations such as, e.g. evolution of a small (sub)system 
interacting with a big system in an equilibrium state (thermostat). For a many- 
particle system, such as considered in the previous section, it is more convenient 
to deal with the obtained in Sec. 2 nonlinear generalized master equation l|19fl 
(rather with its simplified form iJjUJl). For example, as it was shown above, 
equation (|32[1 allows for obtaining the nonlinear equation for a one-particle dis- 
tribution function in the space inhomogeneous case which is not very easy for 
consideration by means of the BBGKY hierarchy pQ . 

Now, we are going to apply the approach of E3] to equation lp?2*|l . Our 
goal is to obtain from (|32[1 the homogeneous equation containing the irrelevant 
part of a distribution function (statistical) operator fi(to) in the "mass" (su- 
per) operator acting on the relevant part of a distribution function (statistical 
operator) f r (t). Such an equation will allow for treating the initial correlations 
(contained in /j(io)) on the equal footing with all other correlations (collisions) 
on any timescale including the initial stage of evolution when initial correla- 
tions matter for sure. It will be the case, because the perturbational expansion 
of a "mass" (super)operator does not lead to the appearance of the "secular 
terms" , proportional to time and violating the perturbational expansion at the 
large (kinetic) timescale. It may be also expected that the divergencies appear- 
ing in the higher order terms in the gas density expansion [21 E] will be more 
conveniently tackled by the "mass" (super)operator expansion. In contrast to 
the HGMEs obtained in j5J j the equation derived below is a nonlinear one 
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and also applicable to the non-conservative systems (with the time-dependent 
Hamiltonians). 

Following 9 , let us identically present the irrelevant part of a distribution 
function (statistical operator) as 

fi(to) = F(t )-f r (t )=Q(t )F(t ) 

= [Q(t )F(t )] F-^top-^t, t )[P(t) + Q(t)}U(t, t )F(t ) 
= C U- 1 (t,t )[f r (t) + f i (t)], 
C a = [QitoMto^F-Hto^U- 1 ^) =T_eiq>[-fd S L(s)}, (54) 

to 

where t/ _1 (i, to) is the backward-in-time evolution operator for the density ma- 
trix F(t), U^ 1 ^, to)U(t, to) = 1 (compare with T_ is the antichronological 
time-ordering operator arranging the time-dependent operators L(s) in such a 
way that time arguments increase from left to right, F (to) is invert to F(to), 
F- 1 (t )F(t Q ) = 1, P(t) + Q(t) = 1. Thus, the additional identity JSU has 
been obtained by multiplying the irrelevant part by unity F~ 1 (to)F(to) (which 
implies the existence of F _1 (io)) and inserting the unities C/ _1 (<, to)U(t, to) = 1 
and P(t) + Q(t) = 1. Therefore, neither divergency (caused by possible vanish- 
ing of F(to)) nor indetermination of the 0/0 type (behaviors of the numerator 
and denominator in F(to)/F(to) = 1 are similar ) may happen. This holds over 
all further (identical) manipulations (see below). 

In l|54|) the following parameter of initial correlations is introduced 

Co = [Qi^Fito^F-^to) = fi(to)[fr(to) + MtoT 1 
= Mto)f- 1 (to)[l + Mt )f- 1 (t )}- 1 

= (l-C )fi(t )f- 1 (t ). (55) 

It is important to note that in (|54|l and i|55|) operator <3(^o) acts only on F(to), 
which is reflected in comprising Q(t )F(t ) into the parenthesis. It follows from 
the fact that f r (t) and fi(to) are the basic quantities we are dealing with in equa- 
tion l|32l) . All functions of dynamical variables, the average values of which we 
can calculate with the help of f r (t) (f r ed{t)) by multiplying equation lIML'l) with 
the corresponding functions (operators) from the right and calculating an aver- 
age value (a trace) , are dependent only on the variables which are not integrated 
off by D (D integrates off all excessive variables in F(t)). Therefore, if we repre- 
sent f r (t) and fi(t ) in ^ as f r (t) = P(t)F(t) and (t ) = Q(t )F(t ), corre- 
spondingly, then the projection operators P(t) = C(t)D and Q(t) = 1 — P(t) in 
these expressions act only on F(t) but not on the functions (if any) to the right 
of them. This is the essence of the reduced description method, when, in order 
to calculate the average values of the functions dependent on a much smaller 
number of variables than the whole distribution function F(t), we actually need 
only the reduced distribution function (density matrix) f re d(t)- 

As it is seen from il53)l . the correlation parameter is a series in /i(to)/ r ^ 1 (*o) 
and, therefore, one may only need a formal existence of the function (operator) 
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/~ 1 (to), which is invert to the relevant distribution function (statistical opera- 
tor) chosen with the help of the appropriate operator P(to). It seems plausible 
that the invert relevant part of a distribution function (statistical operator) de- 
fined in a pointed above sense (uncorrelated part) may always be constructed 
(see E3 an d the next section) . 

Now, we have two equations, i(TT)|) and relating fi(t) to fi(t ). Using 

relations (|30() . finding fi(to) from the pointed equations as a function of f r {t) 
and inserting it into ((32(1 . we obtain the following equation 



- DL°(t)f r (t)+DL'(t)R(t,t )f r (t) 

dP(t) 



l t l I 

DL (t)R(t,t )Jdt U(t,t 

to 



Q(t )L{t ) 



or 



fr(t), (56) 



where the operator R(t,to) is defined as 
R(t,t ) =1 + C(t,t ), 

C(t,t ) = U(t,t ) [1 - C'o^to)}- 1 C U- l {tM), 
C (t,t ) = CqC/^MoMMo), (57) 

Df r {t) = fred{t) and relation (|2"7|l holds. 

Equation is the central result of the paper. We have derived the de- 
sirable homogeneous generalized evolution equation for the relevant part of a 
distribution function (statistical operator). This equation differs from the lin- 
ear homogeneous time-convolution generalized master equation (TC-HGME) 
obtained in [H] from the conventional Nakajima-Zwanzig TC-GME with the 
help of time-independent projection operator technique. Equation 1(56(1 is gen- 
erally a nonlinear one, valid for the systems with time-dependent Hamiltonian 
and holds in both the classical and quantum physics cases if the proper redef- 
inition of the symbols is done and all (super)operators exist (we will address 
the latter problem below). We have not removed any information while deriv- 
ing equation l(56|) . and, therefore, it is exact integra-differential equation which 
accounts for initial correlations and their dynamics through the modification 
of the (super)operator (memory kernel) of ((32(1 acting on the relevant part of 
a distribution function (statistical operator) f r (t). The obtained exact kernel 
of 1(56(1 may serve as a starting point for consecutive perturbation expansions. 
In many cases such expansions of the homogeneous equations (like JHHJ) have 
much broader range of validity than that of the inhomogeneous (like (J22J)) ones, 
when the expansions of the functions (/ r , fi), rather than equation, are involved 
(see also pQ). Equation 1)56(1 reduces to the linear TC-HGME obtained in [H] 
if the (projection) operator P and the Hamiltonian H do not depend on time 
and relations 127(1 . 1(28(1 hold. Equation 1(56(1 is suitable for both, the many-body 
systems, evolution of which is described, naturally, by the nonlinear (like the 
Boltzmann) equations, and the open systems like the ones interacting with the 
thermostat (the linear version of 1(56(1'). How the linear TC-HGME works in the 
classical and quantum physics cases has been demonstrated in and |1()( . 
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The problem of the existence (convergency) of R(t,to) may be raised. The 
function R(t,to) behaves properly at all times. Moreover, the expansion of the 
kernel of (|56|) may result in canceling the pole in the function R(t, to). For linear 
TC-HGME it has been shown in (UJ ^3 m the linear approximation in the small 
density for a space homogeneous dilute gas of classical and quantum particles. 
In such cases there is no problem with the existence of R(t, to). Looking at l|55|l 
(from which the relation (1 — Co) _1 Co = fi(to)/fr(to) follows), one may expect 
that the same will be the case for R(t,to) l|57|l (see below). 

5 Evolution equation accounting for initial cor- 
relations for a dilute gas of particles 

Let us now apply equation i(56|) to a system of N (N ^> 1) identical interacting 
classical particles described by the Liouvillian l(33|) and considered in Sec. 3. 
In order to get an equation for a one-particle distribution function F x (xi , t) 
(i = 1,...,N), we will use the same operator P(t) ((36|) . Thus, taking into 
account (|34|) . (|42|l and 1(43(1 . we obtain from (|56() the following equation for a 
one-particle distribution function 

dF x {x x ,t) 



Of 



= L° 1 Fi(x 1 ,t)+n I dx 2 L 12 F 1 (x 2 ,t)F 1 (x u t) + DL ' C(t,t )f r (t) 

fr(t)- 



DL Jdt R(t,t )U(t,t 

to 



{I?-Ll) + Q{t)L - d -^- 



(58) 



Again, let us consider equation l|58|) in the linear approximation in the den- 
sity parameter 1(46(1 . Given that action of P(t) or D results in the expressions 
proportional to at least the first power in the density n and that all terms on the 
rhs of 1(58(1 (except the first flow term) are already proportional to n (or D), we 
may in the first approximation in n disregard all terms with P(t) in 1)58(1 includ- 
ing such terms in U(t,t ). Then, in this approximation U(t,t ) = U(t,t) and 
is given by (|4*7I) . Therefore, C/ — 1 (t, t )U(t, t n ) in (|57jl may be approximated by 
unity and Co(t, to) by Co- Thus, in the adopted approximation the correlation 
parameter 

C(t,t ) = U(t,t )[f i (t )f- 1 (t )}U- 1 (t,t ) (59) 

as it follows from 1(55(1 and 157(1 . 

To calculate the contribution of initial correlations to the evolution equation 
jSg} , given by (J55J), the definitions for f r (t a ) and fi(t ) iglj should be used. 

Now, proceeding as in Sec. 3., i.e. using (|3~4~|) . i|4*2"|) . l|4^|) and the fact (pointed 
above) that each additional integration over the phase space leads to the addi- 
tional power of n, equation 1(58(1 may be presented in the first approximation in 
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n as 



at 



LlF^xu^ + n I dx 2 L l2 [l + C 12 (t-to)}F 1 (x 2 ,t)F 1 {x 1 ,t) 

t-t 

I 





i / dx 2 L 12 [l + C 12 (t - t )} f°dt 1 e L ^[(L 2 +L' 12 )F 1 (x 2 ,t-t 1 )F 1 (x 1 ,t-t 1 ) 
J o 



where 



dF 1 (x 2 ,t- h) 

H 777 Fi(xi,t - ii)], (60) 

Ot\ 

CWf - in) = e £ "C«-*o) ga^i'^) -z, 13 (t-«o) (g-n 

L X2 = L° 2 +L' 12 , L\ 2 = L\ + Ll and F^Xi) = F 1 (x i ,t ). 

Homogeneous equation I|tj0|) is the central result of this section. It exactly 
describes (in the linear approximation in n) the evolution of a one-particle dis- 
tribution function of a nonideal classical gas of particles with arbitrary space 
inhomogeneity and at any timescale. This equation exactly (in the adopted 
approximation in n) accounts for initial correlations and treats them on the 
equal footing with the collisions between particles and other processes. Initial 
correlations enter the memory kernel (mass operator) governing the evolution 
of the one-particle (reduced) distribution function and do not result in any un- 
desirable inhomogeneous term in the evolution equation l|60|) (compare with 
(O). Contribut ion of initial correlations to the evolution process is given by 
the time-dependent correlation parameter H61|) which describes the evolution 
in time of the normalized two-particle correlation function ^ 9 ?^ x \^f > \ ■ The 
time-evolution of all terms in H60|l is defined by an exact two-particle propa- 
gator G\ 2 (t) = exp(Li 2 t). Again, it is worth stressing that equation (j60|l is a 
nonlinear one, i.e. contains the products of the distribution functions of the 
interacting particles at the same time whereas equation, obtained in |5] from 
the linear TC-HGME in the same, linear in n approximation, is in fact a linear 
equation for F\(x\,t) containing the products of Fi(x\,t) and Fi{x 2 ,to). For 
converting the latter equation into a nonlinear one a coarse-graining over the 
interval t — io <C t re i was needed (on this time interval there is no difference 
between Fi(x 2 ,to) and Fi(x 2 ,t)). 

Let us consider the rhs terms of (|60H and their evolution in time in more de- 
tail. The first term is a conventional flow term caused by the gas inhomogeneity. 
The second term represents the Vlasov self-consistent field modified by initial 
correlations. If C% 2 (t — t ) vanishes (under the action of G\ 2 {t — to)) with time 
at t — to > t cor , we get the conventional Vlasov term. The third term in the rhs 
of (|60|l describes the collisions between particles and the influence of the change 
of the one-particle distribution function on the microscopic scales in space and 
time on the collisions. All these processes are modified by initial correlations 
given by C\ 2 (t — to). Thus, equation (|f>()|) accounts for influence of all pair col- 
lisions and correlations on the dissipative and nondissipative characteristics of 
the nonideal space inhomogeneous gas of particles. 
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Noting that the term with L 2 is of the order of ro jr^ and the term with the 
time-derivative of Fi(x 2 ,t) is ~ t cor /t re i (~ 7), these terms may be disregarded 
(if one is not interested in the nondissipative corrections to the nonideal gas 
dynamics) in the adopted (linear in n) approximation. Then, equation l|6U[l 
takes the simpler form 



mxxut) 

dt 



= L\Fi{xx,t) + n J dx 2 L 12 [l + C 12 (t - t )]F 1 (x 2 ,t)F 1 (x 1 ,t) 
+ n dx 2 L 12 [l + C 12 {t - t )} 



t-t 

I 





x f dtie Ll2tl i 12 F 1 (a; 2 ,i-ti)F 1 ( a; i,i-ti), (62) 



where in the last (collision) term the distribution functions of the collided par- 
ticles may be taken at the same space point ri = r 2 with the adopted accuracy 
~ fo/fh 1 (in the case of plasma, it is not good to do that in the second term 
due to the long-range Coulomb interaction). 

If we consider the kinetic timescale, t — to> t re i 3> t cor , and the two-particle 
correlation function C\ 2 {t— to) (|61|l vanishes at this timescale, equation (|62|l may 
be (for the reasons outlined in the Sec. 3) rewritten in the form of irreversible 
Markovian kinetic equation 



dt 



(>h[ - r ' J:) = L° 1 F 1 (x u t) + n J dx 2 L 12 F 1 (x 2 ,t)F 1 (x u t) 

dx 2 L 12 J dt^ 12 * 1 L 12 F 1 (x 2 ,t)F 1 (x 1 ,t). (63) 




For an inter-particle interaction with the small parameter e <C 1, the two- 
particle propagator exp(Li 2 t) in the collision term of (|63|l may be substituted 
with a "free" propagator exp(L5 2 i) (see (JSOJ) an d equation reduces in the 
second approximation in e to the Vlasov-Landau equation (see, e.g. |14|b In the 
space-homogeneous case equation l|63|l coincides with the Boltzmann equation 
for the momentum distribution function Fi(pi,t). 

Note, that all stages of the evolution of the system under consideration 
may be followed (in the adopted approximation) with the help of equation l|6U|l 
beginning from the initial reversible regime to ^ t ^ t cor . If the reversible terms 

in (|6()J) . related to initial correlations C\ 2 {t — to), vanishes with time and the 

00 

integral J exists, then the evolution of the system may switch from reversible 


to irreversible stage described by, e.g., kinetic equation (|63|l . Thus, we have 
developed an approach to studying the system's evolution which includes the 
evolution in time of initial correlations and does not rely on the RPA or the 
principle of weakening of initial correlations. 
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6 Conclusion 



In this paper, the nonlinear generalized master equations describing the evolu- 
tion of the relevant part of a distribution function (statistical operator) have 
been derived. This new approach is based on the use of the nonlinear time- 
dependent operator P(t) defining the relevant part of a distribution function 
(statistical operator). Generally, this operator is not a projection operator. The 
inhomogeneous nonlinear generalized master equation (|19f) (and may be 

viewed as a generalization of the Nakajima-Zwanzig time-convolution general- 
ized master equation to the case of the time-dependent operators P and Q. For 
the case of time-independent projection operators P and Q this equation reduces 
to the conventional TC-GME. The obtained inhomogeneous nonlinear GME is 
equally useful for deriving both the nonlinear and linear evolution equations for 
the reduced distribution function (statistical operator) of interest in contrast to 
the linear TC-GME which is naturally more convenient for deriving the linear 
evolution equations, e.g. for a subsystem interacting with a thermostat. Using 
the obtained equation, we have derived the inhomogeneous nonlinear equation 
<|49l) for a one-particle distribution function in the linear approximation in the 
small density parameter of a gas of classical particles. This equation, valid for 
arbitrary space inhomogeneity, contains the space and time changes of the one- 
particle distribution function on the microscopic timescale (contributing to the 
nondissipative characteristics of a nonideal gas) and an irrelevant part (a source) 
given by the initial two-particle correlation function. Although the obtained in- 
homogeneous nonlinear GME allows for deriving the closed nonlinear equation 
for a one particle distribution function, the additional assumption like RPA or 
the Bogoliubov principle of weakening of initial correlations is needed in order 
to get either the Vlasov-Landau equation or the Boltzmann kinetic equation. 

For including the initial correlations into consideration we have applied the 
method suggested in and |1U| to converting the inhomogeneous nonlinear 
generalized master equation into the homogeneous form. The obtained ex- 
act homogeneous nonlinear generalized master equation (|56|) for the relevant 
part of a distribution function (statistical operator) describes all stages of the 
(sub)system of interest evolution including the initial stage when the initial cor- 
relations matter. For deriving this equation, which may be used for getting 
both the nonlinear and linear evolution equations for the reduced distribution 
functions, no approximation like RPA or the principle of weakening of initial 
correlations has been used. To test this equation, we have applied it to a space 
inhomogeneous dilute gas of classical particles and obtained in the linear in the 
density approximation a new homogeneous nonlinear evolution equation IjtjUD for 
a one-particle distribution function retaining initial correlations in the memory 
kernel. This equation differs from the one obtained (in the same approxima- 
tion for the space-homogeneous case) in 9 , which in fact is a linear equation 
derived from the linear homogeneous generalized master equation. Equation 
H6Q(I describes all stages of the evolution, accounts for all two-particle corre- 
lations (collisions) and converts to the Boltzmann or Vlasov-Landau equation 
on the appropriate timescale if all initial correlations vanish on this timescale. 
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It is expected, that the next terms in the density parameter expansion (in- 
cluding possible divergencies) will be more easily treated in the framework of 
the obtained homogeneous generalized master equation than by conventional 
approaches [21 IH] - 
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